A new class of exponential propagation iterative methods of Runge-Kutta type (EPIRK)

نویسنده

  • Mayya Tokman
چکیده

We propose a new class of the exponential propagation iterative methods of RungeKutta-type (EPIRK). The EPIRK schemes are exponential integrators that can be competitive with explicit and implicit methods for integration of large stiff systems of ODEs. Introducing the new, more general then previously proposed, ansatz for EPIRK schemes allows for more flexibility in deriving computationally efficient highorder integrators. Recent extension of the theory of B-series to exponential integrators [1] is used to derive classical order conditions for schemes up to order five. An algorithm to systematically solve these conditions is presented and several new fifth order schemes are constructed. Several numerical examples are used to verify the order of the methods is verified and illustrate the performance of the new schemes.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

New Adaptive Exponential Propagation Iterative Methods of Runge-Kutta Type

Exponential integrators have emerged as an efficient alternative to commonly used time-integrators. Recently a new class of exponential propagation iterative methods of Runge-Kutta type (EPIRK) has been introduced [30]. These schemes possess a structure that makes them computationally advantageous compared to other exponential methods. In addition, the general EPIRK formulation offers flexibili...

متن کامل

EPIRK-W and EPIRK-K time discretization methods

Exponential integrators are special time discretization methods where the traditional linear system solves used by implicit schemes are replaced with computing the action of matrix exponential-like functions on a vector. A very general formulation of exponential integrators is offered by the Exponential Propagation Iterative methods of Runge-Kutta type (EPIRK) family of schemes. The use of Jaco...

متن کامل

Solving Parameter Estimation Problems with Discrete Adjoint Exponential Integrators

The solution of inverse problems in a variational setting finds best estimates of the model parameters by minimizing a cost function that penalizes the mismatch between model outputs and observations. The gradients required by the numerical optimization process are computed using adjoint models. Exponential integrators are a promising family of time discretizations for evolutionary partial diff...

متن کامل

Exponential Runge-Kutta Methods for Stiff Kinetic Equations

Abstract. We introduce a class of exponential Runge-Kutta integration methods for kinetic equations. The methods are based on a decomposition of the collision operator into an equilibrium and a non equilibrium part and are exact for relaxation operators of BGK type. For Boltzmann type kinetic equations they work uniformly for a wide range of relaxation times and avoid the solution of nonlinear ...

متن کامل

Runge-Kutta convolution quadrature for operators arising in wave propagation

An error analysis of Runge-Kutta convolution quadrature is presented for a class of nonsectorial operators whose Laplace transform satisfies, besides the standard assumptions of analyticity in a half-plane Re s > σ0 and a polynomial bound O(s 1) there, the stronger polynomial bound O(s2) in convex sectors of the form | arg s| ≤ π/2 − θ < π/2 for θ > 0. The order of convergence of the Runge-Kutt...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • J. Comput. Physics

دوره 230  شماره 

صفحات  -

تاریخ انتشار 2011